Suppression of oscillations by Levy noise 



A. I. 01emskoi a ' b '**, S. S. Borysov b '*, I. A. Shuda b 

""Institute of Applied Physics, Nat. Acad. Sci. of Ukraine, 
58 Petropavlovskaya St., 40030 Sumy, Ukraine 
b Sumy State University, 2 Rimskii-Korsakov St., 40007 Sumy, Ukraine 



Abstract 

We find analytical solution of pair of stochastic equations with arbitrary 
forces and multiplicative Levy noises in a steady-state nonequilibrium case. 
This solution shows that Levy flights suppress always a quasi-periodical mo- 
tion related to the limit cycle. We prove that such suppression is caused by 
that the Levy variation AL ~ (At) 1 / with the exponent a < 2 is always 
negligible in comparison with the Gaussian variation AW ~ {At) 1 / 2 in the 
At —7- limit. Moreover, this difference is shown to remove the problem of 
the calculus choice because related addition to the physical force is of order 
(Atf/ a < At. 
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1. Introduction 

It is known crucial changing in behavior of the systems that display noise- 
induced [U [2] and recurrence [31 H] phase transitions, stochastic resonance 
[SJ E], noise induced pattern formation [H E], noise induced transport P 12] 
etc. is caused by interplay between noise and non-linearity (see Ref. [TO] , 
for review). Noises of different origin can play a constructive role in dy- 
namical behavior such as hopping between multiple stable attractors [TTJ [12] 
and stabilization of the Lorenz attractor near the threshold of its formation 
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[T31 H3]. This type of behavior is inherent in finite systems where exam- 
ples of substantial alteration under effect of intrinsic noises give epidemics 
[E] [E], predator-prey population dynamics [13 [19], opinion dynamics [20J, 
biochemical clocks [2~T] 122], genetic networks [23] . cyclic trapping reactions 
[24J et cetera. 

Above pointed out phase transitions present the simplest case, when joint 
effect of both noise and non-linearity arrives at non-trivial fixed point appear- 
ance only on the phase-plane of the system states. In this consideration, we 
are interested in studying much more complicated situation, when stochastic 
system may display oscillatory behavior related to the limit cycle appearing 
as a result of the Hopf bifurcation [2S1 26J. It has long been conjectured [2Z] 
that in some situations the influence of noise would be sufficient to produce 
cyclic behavior [28J . Moreover, it has been shown that excitable [29] , bistable 
[30] and close to bifurcations [31] systems display oscillation behavior, whose 
adjacency to ideally periodic signal depends resonantly on the noise inten- 
sity [32J (due to this reason, such oscillations were been called coherence 
resonance [29] or stochastic coherence [10J). 

Characteristic peculiarity of above considerations is that all of them are 
restricted by studying the Gaussian noise effect, while such a noise is a special 
case of the Levy stable process (the principle difference of these noises is 
known [33J to consist in the form of the probability distribution that exhibits 
the asymptotic power-law decay in the latter case and decays exponentially 
in the former one). Nowadays, anomalous diffusion processes associated with 
the Levy stable noise are attracting much attention in a vast variety of fields 
not only of natural sciences (physics, biology, earth science, and so on), but 
of social sciences such as risk management, finance, etc. 

In the context of physics, recent investigation [31] has shown that joint 
effect of both non-linearity and Levy noise may cause the occurrence of gen- 
uine phase transitions which relates to a fixed point on the phase-plane of the 
system states. In this connection, natural question arises: may be displayed 
a self-organized quasi-periodical behavior related to the limit cycle by a sys- 
tem driven by the Levy stable noise? This work is devoted to the answer to 
above question within analytical study of two-dimensional stochastic system. 

The paper is organized as follows. In Section [2j we consider pair of 
stochastic equations with arbitrary forces and multiplicative Levy noises to 
obtain their analytical solution in a steady-state nonequilibrium case. This 
allows us to conclude in Section [3] that opposite to the Gaussian noises the 
Levy flights suppress always a quasi-periodical motion related to the limit 
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cycle. Since equation, governing behavior of stochastic system driven by 
multiplicative Levy stable noise, are very complicated [35J and moreover 
their derivation is now in progress [36], we complete our consideration with 
Appendix A containing details of derivation of the Fokker-Planck equation. 
Moreover, to demonstrate that a closed consideration of the Levy processes is 
achieved only within the Fourier representation we set forth a scheme related 
to the appropriate stochastic space in Appendix B. 



2. Statistical picture of limit cycle 

According to the theorem of central manifold [25], to achieve a closed 
description of a limit cycle it is enough to use only two degrees of freedom 
related to some stochastic variables X iy % = 1, 2. In this way, stochastic evo- 
lution of the system under investigation is defined by the Langevin equations 
[37] 

dXi = f4t + gi dLi, i = l,2 (1) 

with arbitrary forces /j = fi(xi,x 2 ) and noise amplitudes Qi = gi(xi,x 2 ) 
being functions of both variables Xj, i = 1,2; stochastic terms are related to 
the Levy stable processes L; = Lj(i). Within the Ito calculus, these processes 
are determined by the elementary characteristic function 



c 



with increments Li = Li(k\, k 2 ; xi, Xq) whose expression [35J 
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& = ih (fi + 7tfi) - {rriigikifie-W® ^ \m jgj k,fl e' 1 ^ (3) 

3=1 



follows from Eq.(A.25). Hereafter, we use asymmetry angles ipi and moduli 
rrii defined by the equalities 

tan [<Pi(a)} = ^sgn^fcj) tan(7ra/2), 

(4) 



m a % = ^/l + /32tan 2 (7ra/2); 



everywhere, the Levy index a G (0, 2) characterizes the asymptotic tail 
x i ^ a+1 ^ of the Levy stable distribution at 1 ^ a < 2 (the case a = 2 relates to 
the Gaussian distribution), parameters /3j G [—!,+!] define the distribution 
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asymmetry, location parameters — oo < 7$ < +00 denote the mean values of 
stochastic variables Xi at a > 1, and the angular brackets denote averaging 
over Levy noises. 

As is shown in Appendix A, the Fourier transformed probability distri- 
bution function 



+00 



P{k 1 ,k 2 ;t) = T{P{x 1 ,x 2 )}{k 1 ,k 2 ;t) :-- 



is governed by the Fokker-Planck equation 

2 



dxidx 2 P{ Xl ,x 2 ;t)e iiklXl+k2X2) 

(5) 



dP 
~dt 



i=l 



P. (6) 



Characteristically, being Fourier transformed, r.h.s. of this equation depends 
on the wave vector components ki and k 2 , while both forces /j = fi{xi,x 2 ) 
and multiplicative noise amplitudes ^ = 5 f i(xi,x 2 ) are dependent on the 
coordinate components X\ and x 2 . 

According to the continuity equation (A.23), components of the steady- 
state probability flux are obeyed to the condition ^dJj/dxj = which 
means the first component J\ = Ji(x 2 ) is a function of the only variable 
x 2 , and vice- versa for the second component J 2 = J 2 {x\). Then, within the 
Fourier representation, the system behaviour is defined by the equations 



(fi +0i7i) +iK#i|fe ^(^l^lf- 2 ^ 



x ^m^hfl e~ ilp1 ^ + \m 2 g 2 k 2 fl e' iip2 (^ } P = 2vr J 1 (k 2 )5(k 1 



(7) 



{(/2 + £272) + i\m 2 g 2 \*e ^(S) \k 2 \ a 2 A; 2 



x f|mic/ifci| f e" ivi (t) 



+ \m 2 g 2 k 2 



a ■ ( a \ 



P = 2irJ 2 (k 1 )5(k 2 ). 



Since the pair of these equations determines a single distribution function 
P(ki, k 2 ), the consistency condition 



[(/1 + <7i7i) + ie-^^lmi^iri^ir" 2 ^] S{k 2 )J 2 {h) 
[(/ 2 + ^272) + ie-^^^llfer 2 ^] SikJMh) 



(9) 
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should be kept to restrict the choice of the probability flux components J\{k 2 ) 
and J 2 {k\). 

Multiplying Eq.(7) by the factor |m 2 g 2 | ^e~ lv,2 ( 2) and Eq.(8) by \migi\^e~ 1<Pl ^ 
and then subtracting results, one obtains 

F + i|mim 2 ^ 2 |f e _i h(t)+^(f )] 



x 



\mtgikil 2 e 



\m 2 g 2 k 2 \* e~ ilp2 ( 



2tt 



Ji{k 2 )8(ki)\m 2 g 2 \ 



J2(k 1 )5(k 2 )\m 1 g 1 \^e- i ^ 



where one denotes 

F = (/ 1 + 7 1 ^ 1 )|m 2 ^|te- i ^(t; 



- {h + 72^2) |miflfi|ae" 



(10) 
(11) 



The equation (10) yields the explicit form of the probability distribution 
function 

u( Tdfe J 1 (fc 2 )|m 2 ^ 2 |fe- i [ fc ^ + ^(f)] 

-T (£1, X 2 J — / 

—00 
+00 



2tt F 2 


- ibi 


2 m 2 # 2 a e" 


-i V2 (a)\k 2 \ a - 2 k 2 










[fcl!Bl+^l(f)] 


(12) 




^1) 


l^i^il^e" 1 




2vr F a 


+ %2| 


2 mi^i a e _ 


-' 1 ^i a )\ki\ a - 2 kx 





where effective forces F 12 are determined by Eq.( 11 ) at m 21 = 1 and ip 2 ,i = 0. 

In the case of constant values of the probability flux within the state 
space Xi,x 2 , the Fourier transforms related are J^k-i) = 2nji 0) 5(k 2 ) and 
^2(^1) = 2irJ 2 S(ki) with jj® = const. Then, the consistency condition (9) 
takes the form (fi + (7171) J 2 °^ = (f 2 + g 2 "y 2 ) j{ , the effective force (11 ) is 
Fo = {h + 7i0i) Iflbl* ~ (/a + 72^2) |0i|*, and the probability density (|12|) 
reads 



P 



J l 1^2 1 



J 2 I Pi I 



(13) 



To create a limit cycle this distribution function should diverges on a closed 
curve, so that the effective force equals F = 0. Together with the consistency 
condition, this equation gives 



4 0) _ h + 7igi 

Ji 0) /2 + 7202 
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(14) 



But these equalities mean that the numerator of the probability density (13) 
disappears also. As a result, we conclude the limit cycle creation is impossible 
for a stationary non-equilibrium state with both probability flux components 
Ji(xi,x 2 ) and J 2 (xi, x 2 ) being constant. 



To calculate integrals in Eq.(12) for arbitrary dependencies J\{k 2 ) and 
J 2 {ki) it is convenient to write \k\ = sgn(k)k = e me (~ k )k where 9(k) denotes 
the Heaviside step function. Then, one has \k\ a ~ 2 k = e ~ 17r6 ' (_fc )( 2 ~ a )/c Q ~ 1 , and 



the pole points of integrands in Eq.(12) are expressed with the equality 
F h2 \ ^ 

i i a I 

y I ^ I — / 

(15) 



|wil,201,2r|02,l| 



x exp < i 



<p lt 2 (a) + (2 - a)7r6(-^K h2 ) + (7r/2)sgn(^ 1;2 ) 

a — 1 



Due to sign-changing term (7r/2)sgn(Q : i^i i 2) in the exponent the Ki >2 poles 
are located on opposite half-planes of complex variables k\ )2 . Making use of 
the power series expansion 

k a-i = Ra -i 1 1 + ^ZJ}_ j ^ K^ 1 + {a-l)K a ~ 2 {k- K) (16) 



allows us to reduce the integrands in Eq.( 12 ) to a pole form. However, we can 
not close the integration contours around both upper and lower complex half- 
planes of the k variable since integrands related contain absolute magnitudes. 

To find the integrals needed let us specify a contribution that gives the 
pole located on the upper half-plane of the complex number k. With this aim, 
we divide this half-plane into two parts related to the positive and negative 
values of the real part of k. As shows Figure [TJ integrals in Eq.(12) can be 
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rewritten as follows: 



+00 



f(k) 
k-K 



dk 



ABC 



f(k) 
k-K 

f(k) 
k-K 



DEF 



f(k) 
k-K 



dk + 



k-K 



dk 



Wc 



FD 



k-K 



dk + 



k-K' 



dk 



AB 



DE 



dk- 



dk 



WC 



WF 



f(k) 
k-K 

f(k) 
k-K 

f(k) 
k-K 



dk + 



CA 



dk + 



f(k) 
k-K 

f(k) 
k-K 



dk 



dk 



FD 



dk + 



f{k) 
k-K 



dk 



ABC 



DEF 



f{k) 
k-K 



dk + 



f(k) 
k-K 



dk 



\CA 



EF 



(17) 



With tending radiuses of the arcs CA and EF to infinity, both integrals in 
the last square brackets disappear. On the other hand, when both half-axes 
BC and FD tend one to another, one has J BC = — j FD , so that terms in the 
square brackets standing before are cancelled also. Moreover, the integral 
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over the contour DEF equals zero because this contour does not envelop 
any pole. As a result, we obtain 



+00 



f(k) 
k-K 



dk 



k-K 



dk = sgn(%K)2nif(K) 



(18) 



ABC 



where the last equality is due to the residue theorem. 



Finally, making use of the Cauchy integral (|18|) yields the probability 

(19) 

J21 (K h2 ) 



distribution (12) in the form 



where one denotes 

Pl,2 = 



a a(3 — a) 5 

(a - 1)1^2,1 1 2(a - 1J 1^1,201,2 1 2(a " 1) 

j .2 == 3 —^'Pi,2 (%) + [sgn (3^ 1)2 ) + 26 (-^1,2)] • 



(20) 



a — 1 



2 a - 1 



3. Discussion 



Analytical consideration developed in previous Section allowed us to ob- 



tain the probability distribution function (19) that describes behaviour of 



nonequilibrium steady-state stochastic system driven by the Levy multiplica- 
tive noise with two degrees of freedom. Recently, we have studied conditions 
of the limit cycle creation in stochastic Lorenz-type systems driven by Gaus- 
sian noises [38J. Noise induced resonance was found analytically to appear 
in non-equilibrium steady state if the fastest variations displays a principle 
variable which is coupled with two different degrees of freedom or more. The 
condition of this resonance appearance is expressed formally in divergence of 
the probability distribution function being inverse proportional to an effec- 
tive force type of (11) - when this force vanishes on a closed curve of phase 



plane, the system evolves along this cycle with diverging probability density. 



In opposite to such a dependence, the distribution function (19) contains 



the effective force ( 11 ) in positive power (2 — a) /(a — 1) only. To this end, we 



can conclude the Levy flights suppress always a quasi-periodical motion re- 
lated to the limit cycle. That is main result of our consideration. The corner 
stone of the difference between stochastic systems driven by the Levy and 
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Gaussian noises is that the Levy variation AL ~ (At) 1 /" with the exponent 
a < 2 is negligible in comparison with the Gaussian variation AW ~ (At) 1 / 2 
in the At — >■ limit. 

It is interesting to note that above difference removes the problem of the 
calculus choice [H [37] . This problem is known to be caused by irregularity of 
the time dependence X(t) of stochastic variable (for the sake of simplicity, 
we consider one- dimensional case again). Hence, in the integral of equation 
of motion ((!]) 

t L(t) 

X(t) = j f(x(t'))dt' + j g(x(t'))dL(t') (21) 

L(0) 

we should take the noise amplitude g(x(t')) at the time moment 

t' = t' + XAt'; X e [0, 1], At' (22) 

that does not coincide with the integration time t! due to a parameter A G 
[0, 1] whose value fixes calculus choice (for example, the magnitude A = 1/2 



relates to the Stratonovich case) [HE?]. With accounting equations (22) and 
([I]), one obtains 

g(x(t)) ~ g(x(t)) + Xg'(x(t))AX(t) 

~ g(x(t)) + Xg'(x(t))f(x(t))At + Xg' (x(t))g(x(t))AL(t) 

where primes denote differentiation over argument x. Being inserted into 



Eq.(21), the first term in the last line of Eq.(23) relates to usual case of the 
Ito calculus. Corresponding insertion of the second term gives an addition 
whose order AL ■ At ~ (At) 1+ ^/ a ^ At is higher than one for the previous 
term (such a situation is inherent in the Gaussian case as well). Finally, after 



insertion of the last term of Eq.(23) the last integrand in Eq.(21) obtains an 



addition of order (AL) ~ (At) ' a . In special case of the Gaussian noise 
(a = 2), the order 2/ a of this addition coincides with the same in the first 



integrand of Eq.(21), that is resulted in addition Xg(x)g'(x) to the physical 
force f(x). Principally different situation is realized for the Levy stable 
process, when the index a < 2 and above addition should be suppressed in 
comparison with physical force because (At) 2//a <C At. 



9 



Appendix A. Derivation of Fokker-Planck equation for the Levy 
multiplicative noises 



Following to the line of Ref. [35J, we start with consideration of one- 
dimensional Levy process X(t) whose Chapman-Kolmogorov equation 

p(x,t + dt\x ,t ) = J dy p(x,t + dt\y,t)p(y,t\x ,t ) (A.l) 

connects transition probabilities taken in intermediate positions y related to 
time t. According to the definition 

p(k, t + dt\y, t) := e d ^( fc ' d *l^), (A.2) 

the inverse Fourier transform 

/dk 
—e- ik (*-y)p(h,t + dt\y,t) (A.3) 

is expressed in terms of the elementary cumulant dKx(k, dt\y, t) of the char- 
acteristic function of the stochastic process X(t). Then, with using the 
dt — > limit and the identity 

p(x,t\x ,t ) = I dy p{y,t\x ,to) j ^~' Mx ~ y \ (A.4) 



the equation (A.l) arrives at the chain of equalities: 

p(x, t + dt\x , t ) - p(x, t\x , t ) 
= J dy pfatM J ^ e - ife M[e d ^^l^-l] 

~ J dy p(y,t\x ,t ) J ^e- ik ^dK x (k,dt\y,t) (A ' 5) 



dy dK x (x - y,t)p(y,t\x ,t ) = dK x (x,t) *p(x,t\x ,t ). 
Here, * denotes the convolution of the inverse Fourier transform 

/dk 
—dK x (k,dt\y,t)e- ik ^-y\ (A.6) 

As a result, with accounting the definition 

dK x (x,t) := C(x)dt, (A.7) 
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the equalities (A. 5) yield the Fokker-Planck equation 
d 

— p(x, t\x , t ) = £(x) *p(x, t\x , t ). (A.8) 

To obtain the explicit form of the increment C(x) let us consider initially 
the Levy process L(t) itself. The elementary characteristic function related 

is determined by the cumulant 

dK L (k,dt\y,t) := A(k)dt (A.10) 

with the Levy increment [31] 

A(k) = ik-f - D\mk\ a e-'^ {a) (A.ll) 

where asymmetry angle cp and modulus m are determined by Eqs. Q. The 
elementary characteristic function of the principle process dX = fdt + gdL 
is written as follows: 



c 



dK x (k,dt\y,t) ._ /tfkdX^ _ Q ikfdt ^(fc^dZ^ _ e ikfdt e dK L (gk,dt\y,t) ^ 12 ) 



where Eq.(A.9) is taken into account. Similarly to the definition (A. 10), the 



elementary cumulant 

dK x (k, dt\y, t) := C{k, x)dt (A.13) 

is determined by the increment 

C(k,x) = ikf(x) + A(g(x)k) (A. 14) 

whose explicit form reads [2U EE] 

£{k, x) = ik [f{x) + ig {x)\ - \mg(x)k\ a e~' l ^ a) . (A.15) 

Hereafter, we renormalize the noise amplitude g(x) to suppress the scale 
factor D. 

The probability distribution function 

P(x,t) = J dx p(x,t\xo,to)P(xo,t ) (A. 16) 

11 



is determined by the Fokker-Planck equation 

;P(x, t) = C(x) * P(x, t) = J dy £(x — y, x)P(y, t) 



d_ 
dt 



(A.17) 



following from Eq.(A.8). After using the Fourier transform 



P(k, t) = F{P(y, t)}(k, t)= dy P{y, t)e 



^iky 



this equation takes the convenient form 





Of 



P(k,t) = C(k,x)P(k,t) 



(A. 18) 



(A.19) 



with the kernel (A.15). It is worthwhile to note this kernel, being the Fourier 



transform inverse to Eq.(A.6) with respect to the coordinate difference x — y, 



depends on the coordinate x through both force f(x) and multiplicative noise 
amplitude g(x). 

With accounting the relation 



0° 



d\x\ 



-h{x) 



-i 



k\ a h(k) 



(A.20) 



for the Riesz derivative with respect to an arbitrary function h(x), Eqs. 



(A.19) and (A.15) arrive at the following form of the fractional Fokker-Planck 



equation [351 ES 



Of 



P(x,t) 



d_ 

dx 



[f{x) + 19 (x)]P{x,t) 



+ 



aa— 1 



\g{x)\ a P(x,t). 



(A.21) 



In symbolic form, many-dimensional generalization of this equation for a 
symmetric Levy flight reads: 



d 



Of 



P(x, t) 



-A : g(x)g(x)] a/2 P(x,t). 



-V[f(x)+g(x)-^P(x,t)- i 

(A.22) 

Here, every dot denotes the summation over indexes % = 1, 2 and the axes 
X\, X2 forming pseudovector x are chosen so that the noise amplitude matrix 



12 



g takes the diagonal form = g^S^ whose elements form the pseudovector 
g. In the component form, one has the continuity equation 



at 

with the probability flux 



8 



x) 







(A.23) 



dxl 







dx« 



[ gi {x) gj {x)]i } P{x) 



(A.24) 

In generalized case of non-symmetric Levy flights, the Fourier transforms of 
the flux components are written in the explicit form 



Jl 



(/i + 0i7i)+i|"H0i| ? e ^i^lhf 2 h 
mxgxhfl e~ ivn (^) + \m 2 g 2 k 2 \^ e~^ 2 ^ 



P. 



P 



(A.25) 



h = { {f'2 + 9212) + i\m 2 g 2 \ze M 2) |fc 2 | 2 2 k 2 
Im^kil^ e~ iipi ^"> + \m 2 g 2 k 2 \^ e~ iip2 (^ 

where we use the asymmetry parameters Q. 



Appendix B. Consideration of the Levy processes within direct 
stochastic space 

After inverse Fourier transformation, the components ([7]) and ^ of the 
stationary probability flux are written as follows: 



Cfi +9ili) + 
(/a + #272) + 



df- 1 






S— 1 

dxl 










(_B 


—— 1 
9x 2 2 




\ dx 



d 



(9192) 



Mi) 2 + 







dx- 



92 



p = 4°\x 2 ), 

P = J { 2 \x 1 ). 



13 



,2 d?- 1 



Acting by the g 2 2 Q „ 1 operator on the first of these equations and the 



® x 2 
a a — l 

gl d \_ x operator - on the second, one obtains 



9i 



ch 1 82 1 
d%2 8x1 



a Q2 1 8 2 ^ 



01 



dxl ® x i 



8 

8x\ 



92 



8 



9i 



8 



8Xr 



(9l92) 



P 



n. 



— — 1 
dx^ 



(929i) 2 + 



4 0) - (A + <7i7i) -P 
9 



(B.l) 



0i 



5i 



J. 



(0) 



<9x 2 ' 2 

(A + 0272) p 



Subtracting above equalities term by term, one arrives at the fractional dif- 
ferential equation 



(A + gill) 92 



a 8 



--1 



-1 



(A + 0272) g'i 



a 8 2 1 



-1 



P + G(xi,x 2 )P 



92 



wzrJr (x 2 ) - gl 

n - 1 - 



(B.2) 



4 0) (xi) 



a_i ^2 

dxi 



where one denotes the function 



G(x!,x 2 ) = gi 



8x2 8x1 



8 



a Ql 



-gi 



dxl & X 2 



dx\ 
8 

dxi 



9i 



(929i) 



8 



dx- 



(9192) 



8 



dx- 



92 



(B.3) 



+ 



9i 



8x2 



:i (A + gxii) - 9i — «zi (A + 0272) 

dxl 



For the Gauss processes (a = 2), the differential equation (B.2) is reduced to 
the algebraic one to give the probability distribution function that has been 
found in our previous work [38]. However, in general case a < 2, solution of 
the fractional differential equation (B.2) arrives at a complicated problem, 
so that we are obliged to use the Fourier representation in Section [2] 
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It is worthwhile to note finally the consistency condition (|9| takes the 
form 



d 

dxi 
d 



(fi + giii) - 

{f 2 + 92l2) - 



d 



dxi 
d 
dxo 



\9i\ 



02 



4°\xi) 



(B.4) 



within the inverse Fourier representation where one takes tpi = and m$ = 1, 
for the simplicity. The equation (B.4) connects explicitly the probability flux 
components J2, 1(^1,2), being arbitrary functions, with given dependencies 
of both forces /i(xi,x 2 ), f^{x\,X2) and multiplicative amplitudes gi(xi,x 2 ), 
92(^1, x 2 ), respectively. 
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